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ABSTRACT 

Context. Spectropolarimetry of AGNs is a powerful tool for studying the structure and kinematics of the inner regions of quasars. 
Aims. We wish to investigate the effects of various AGN scattering region geometries on the polarized flux. 

Methods. We introduce a new, publicly available Monte Carlo radiative transfer code. Stokes, which models polarization induced by 
scattering off free electrons and dust grains. We model a variety of regions in AGNs. 

Results. We find that the shape of the funnel of the dusty torus has a significant impact on the polarization efficiency. A compact 
torus with a steep inner surface scatters more light toward type-2 viewing angles than a large torus of the same half-opening angle, 
Oq. For 6*0 < 53°, the scattered light is polarized perpendicularly to the symmetry axis, whilst for So > 60° it is polarized parallel 
to the symmetry axis. In between these intervals the orientation of the polarization depends on the viewing angle. The degree of 
polarization ranges between 0% and 20% and is wavelength independent for a large range of Oq. Observed wavelength-independent 
optical and near-UV polarization thus does not necessarily imply electron scattering. Spectropolarimetry at rest-frame wavelengths 
less than 2500 A may distinguish between dust and electron scattering but is not conclusive in all cases. For polar dust, scattering 
spectra are reddened for type-1 viewing angles, and made bluer for type-2 viewing angles. Polar electron-scattering cones are very 
efficient polarizers at type-2 viewing angles, whilst the polarized flux of the torus is weak. 

Conclusions. We predict that the net polarization of Seyfert-2 galaxies decreases with luminosity, and conclude that the degree of 
polarization should be correlated with the relative strength of the thermal IR flux. We find that a flattened, equatorial, electron- 
scattering disk, of relatively low optical depth, reproduces type-1 polarization. This is insensitive to the exact geometry, but the 
observed polarization requires a limited range of optical depth. 

Key words. Galaxies: active - Polarization - Radiative transfer ~ Scattering - Dust 



^ 1. Introduction entati on effects (Lawrence & Elvis 119821 De Zotti & Gaskell 

119851 ). Since then, the dusty-torus model has become the stan- 

, One of the foremost problems in AGN research is that the in- jard unified model (see Antonucci 1993^ dividing AGNs into 

nermost regions of AGNs cannot be resolved in the optical and sub-types: "type-1" AGNs which are seen close to face-on, 

UV with current technology. However, the light of AGNs is po- ^nd "type-2" AGNs which are seen close to edge-on. In type-1 

larized over a broad wavelength range, and this allows us to put ^GNs the central energy source and its surroundings (e.g., the 

important constraints on the geometry of the emitting and scat- blr) ^an be seen, whilst in type-2 AGNs the torus blocks our 

tering regions. Spectropolarimetric observations giving the de- direct view of these inner regions. While this obscuration and 

tailed wavelength dependence of the polarized flux give further ^^e IR emission from the torus ai-e the most obvious eff^ects of 

clues to the nature of the polarizing mechanism. the torus, scattering from the dust will add polarized flux. The 

Our inferences of the innermost sti-uctures of AG Ns have so polarization spectrum of an optically thick dusty torus has been 

far been obtained indirectly. Rowan-Robinson ( I1977D suggested the subj ect of several modeling projects (Kartje [19951 Wolf & 

that AGNs are surrounded by a dusty torus and in the same paper Henning I 19991 Watanabe et al. l2003l) . 
he gives a suggestion by M. V. Penston that Seyfert 2 galaxies 

are seen close to edge-on so that the active nucleus is obscured ^hen Dibai & Shakhovskoy (1966) and Walker (1966) dis- 

by the torus. Supportfor this picture came from the important ^^^^^^^^ optj^^l polarization of AGNs, it was initially taken to 

discovery by Keel ([1983 that Seyfert 1 galaxies (active galax- ^e evidence of optical synchrotron emission, since synchrotron 

les showing a broad-hne region; BLR) are preferentially seen radiation has a high intrinsic polarization. However, Angel et 

face-on. Keel 0980) also investigated absorption effects inside ^975, fo^nj the Balmer lines in NGC 1068 to be polarized 

the host galaxies and emphasized the need of additional nuclear similarly to the continuum, thus implying that scattering was re- 

absorption in Seyfert galaxies with respect to normal spirals, sponsible for the polarization of both the lines and continuum. 

Keel's work led to further confirmation of the importance of on- ^he difference they found in polarization between the narrow- 

line region (NLR) and BLR places the scattering region outside 

Send offprint requests to: Rene W. Goosmann the BLR, but inside the NLR 
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When light is scattered, the angle of polarization depends 
on the direction of the last scattering, so one expects the angle of 
polarization to be related to the structure of the AGN. Stockman, 
Angel, & Miley ( 1979 ) made the seminal discovery that for low- 
polarization, high optical luminosity, radio-loud AGNs, the op- 
tical polarization position angles tend to align parallel to the 
large-scale radio structure. Although they interpreted this as a 
consequence of optical synchrotron emission, they also sug- 
gested that polarization from an optically-thin, non-spherically- 
symmetric scattering region near the source of optical radiation 
was another possibility. 

Antonucci ( |1982| l pointed out that whilst many radio galax- 
ies showed a similar parallel alignment of the polarization and 
radio axes, there was, unexpectedly, a population showing a per- 
pendicular relationship. It was subsequently shown (Antonucci 
I1983I I that relatively-radio-quiet Seyfert galaxies show a similar 
dichotomy between the predominantly, but not exclusively, par- 
allel polarization in face-on type-1 Seyferts and the perpendic- 
ular polarization of type-2 Seyferts (see Antonucci 1993, 20021 
for reviews). These discoveries made a synchrotron origin of the 
polarization much less likely. 

Polarization perpendicular to the axis of symmetry is eas- 
ily produced by scattering off material close to the axis. There 
is good observational evidence for the existence of ionization 
cones along the polar axis in numerous objects (see Kinney et 
al. |l991| and reference therein). Polar scattering has been partic- 
ularly well studied in the Seyfert- 1 galaxy NGC 1068. Antonucci 
& Miller (|1985l l made the key discovery that the polarized-flux 
spectrum can offer a periscope view of type-2 AGNs because 
much of the polarized flux originates inside the torus. Detailed 
HST polarimetry has revealed the polarization structure of the 
ionization cones (see Capetti et al. 1995a,b; Kishimoto ,1999) . 

The detection of a hidden broad-line region in NGC 1068 by 
Antonucci & Miller (|19851 l was of great importance for AGN 
research since it provided strong support for the unified theo- 
ries of AGN activity. More hidden type-1 nuclei have subse- 
quently been found by analysis of their polarized-flux spectra 
(see e.g.. Miller & Goodrich, 1990; Tran, Mifler, & Kay, 1992; 
Hines & Wills, 1993; Kay, 1994; Heisler, Lumsden, & Bailey, 
1997; Tran, 2001; Smith et al., 2004). Similar work on the ra- 
dio galaxy 3C 321 was done by Young et al. (1996), and Tran 
et al. ( 1 19991 1 could identify an active nucleus inside an ultra- 
luminous infra-red galaxy using spectropolarimetry. Recently, 
hidden type- 1 nuclei have also been found in five type-2 quasar 
candidates (Zakamska et al. 12005 l l. 

The new generation of large telescopes is delivering spec- 
tropolarimetry of emission line profiles with good velocity reso- 
lution (see, for example, the spectropolarimetry of the Seyfert 
1.5 galaxy NGC 4151 presented by Martel 1998, and the at- 
las of spectropolarimetry of Seyfert galaxies presented by Smith 
et al. 12002 J . Examination of velocity-dependent polarization of 
emission lines promises to reveal valuable information about the 
geometry of the BLR (Smith et al. 120051 1. Similarly, spectropo- 
larimetry of quasar absorption lines helps constrain the geometry 
of broad absorption line QSOs (Goodrich & Mifler 1995, Cohen 
et al. 1995, Hines & WOls 1995, Ogle et al. 1999). 

In order to understand these many facets of AGN polariza- 
tion, and their implications for the underlying geometry, theo- 
retical modeling is necessary. Analytical approaches to radiative 
transfer that have been carried out so far are generally limited to 
the consideration of single-scattering models. Computer simula- 
tions are needed to investigate multiple-scatterings. In this paper 



we describe Stokes, a new general-purpose, publicly-availableQ 
Monte Carlo code for modeling wavelength-dependent polariza- 
tion in a wide variety of scenarios, and we present some results 
of our study of AGN polarization. 

In this paper we confine ourselves to using Stokes for calcu- 
lating the polarization imprints of basic constituents of the uni- 
fied scheme. We compute the polarization spectrum of dusty tori 
with various geometries and opening angles, and we consider 
scattering in polar cones and electron disks. We investigate the 
effects of geometrical shape and optical depth of given regions. 
None of our models are intended to reproduce observational po- 
larimetric data for any specific object. Rather, we want to inves- 
tigate general constraints on the scattering regions and the geom- 
etry of AGNs. We discuss consequences for the observed polar- 
ization dichotomy between type-1 and type-2 objects. We leave 
aside the question of interactions between different types of scat- 
tering regions for paper II (Goosmann & Gaskell, in preparation) 
where we also conduct more detailed modeling of AGNs in the 
unified scheme. 

The present paper is organized as follows: in section 2 we 
summarize modeling of optical and UV polarization of AGN 
and the main results obtained previously. Section 3 describes 
our code Stokes. In section 4 we present modeling results for 
equatorial, toroidal dust distributions. Section 5 is dedicated to 
electron and dust scattering in polar double-cones. In section 6 
we investigate the polarization signature of equatorial regions 
for electron scattering. Our results are discussed in section 7 and 
we give some conclusions in section 8. 

2. Previous codes and modeling 

In this section, we briefly summarize some recent AGN polariza- 
tion modeling codes which we will compare our Stokes model- 
ing with. 

Young et al. (IT9951 [T996] l, Packham (1997), and Young 
(2000) developed an analytical radiative transfer model, the 
Generic Scattering Model (GSM), for polarization modeling. 
The model is based on the unified AGN model. Extended emis- 
sion regions can be defined, and scattering processes as well as 
dichroic absorption are considered. The modeled geometries in- 
clude toroidal, disk-like, and conical regions of dust and free 
electrons. For scattering material in motion, Doppler effects are 
included. The model is fairly effective in reproducing spectropo- 
larimetric data of Seyfert galaxies (see, for example. Young et 
al. 1 1 9991 for Mrk 509). In particular, it reproduces variations of 
the polarization across broad emission lines (Smith et al. 2005Jl. 
The model is semi-analytical and therefore does not take multi- 
ple scattering into account. 

Wolf & Henning ( 1999) present a Monte-Carlo code used to 
compute the polarization obtained by scattering inside axisym- 
metric regions. They consider dust and electron scattering for 
polar double-cones and equatorial tori. In the Monte-Carlo ap- 
proach two or more of such components can be combined and 
the resulting polarization spectra are modeled for various incli- 
nations of the system. Aside from spectropolarimetric model- 
ing, the code by Wolf & Henning ( 1999) can also produce po- 
larization images, which are provided for various torus geome- 
tries. An important element in this code is that multiple scatter- 
ing, which becomes important for optical depths > 0.1, is con- 
sidered accurately by including the dependence of the scatter- 
ing angle and the polarization of a scattered photon on its inci- 
dent Stokes vector. For dust scattering, two different grain size 

' littp://www.stokes-program.info/ 
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distributions were examined: one parameterization representing 
Galactic dust, and the other favoring larger grains. 

Kartje d 19951 1 also developed a Monte Carlo Code and inves- 
tigated quasar schemes with either a torus geometry or conical 
stratified winds along the polar axis. In addition to polarization 
by scattering, he also considers polarization by dichroic extinc- 
tion due to magnetically-aligned dust grains. For a simple uni- 
fied torus model he finds that the dominant parameter of the po- 
larization, P, is the torus half-opening angle: for type-2 objects 
one can find significant polarization (up to 30%) with a posi- 
tion angle directed perpendicular to the axis of symmetry; for 
type-1 objects P is negligible. Kartje obtains an important result 
when he investigates conical stratified-wind regions containing 
free electrons closer to the central source and dust farther out: 
the amount of polarization ranges between 0% and 13%, match- 
ing observed values, and the direction of the ^-vector depends 
on the viewing angle in a manner that agrees with the above- 
mentioned type-l/type-2 dichotomy. The polarization percent- 
age can be increased if there is magnetic alignment of dust 
grains, but the general dependence of P on the viewing angle 
seems to be a geometrical effect. 

Another Monte-Carlo polarization code is presented by 
Watanabe et al. (20031). It is applied to modeling of optical 
and near-infrared spectropolarimetric data of the type-2 Seyfert 
galaxies Mrk 463E, Mrk 1210, NGC 1068, and NGC 4388. The 
code contains electron and dust scattering routines quite similar 
to those used by Wolf & Henning (|1999l l. It considers multiple 
scattering and dichroic absorption in dusty tori, spheres as well 
as electron and dust scattering in double-conical regions. The ab- 
sorption and scattering properties of the dust are carefully calcu- 
lated by Mie theory. Watanabe et al. (2003 ) examine wavelength- 
dependent polarization properties for diff'erent geometries over a 
broad-wavelength range and give constrains about possible scat- 
tering components within the objects they observed. They con- 
clude that a combination of dust and electron scattering in po- 
lar regions can reproduce the optical polarization properties of 
Mrk463EandMrk 1210. The slope of optical polarization NGC 
1068 is almost flat favoring electron scattering as the dominant 
polarizing process. For the near-infrared range polarization of 
these objects can be modeled by dichroic absorption of aligned 
dust grains in a torus. However, scattering off Galactic dust in 
a torus cannot simultaneously reproduce the near-infrared polar- 
ization and the total flux. Watanabe et al. (2003) hence suggest 
that the grain size composition of AGNs might be different from 
our Galaxy. 

This list of previous polarization modeling is not exhaus- 
tive. For example, Blaes & Agol 119961 and Agol & Blaes 



11996 ' have presented modeling of the wavelength-dependent 
polarization signature of accretion disks at the Lyman limit, 
and Kishimoto (1996) modeled polarization due to electron- 
scattering off clumpy media in polar regions of AGN. Also, a 
new Monte-Carlo model, which includes polarization transfer 
for the continuum and for broad quasar absorption lines, was 
recently presented by Wang, Wang, & Wang (2006). We have 
restricted this brief review to recent modeling of polarization by 
dust and electron scattering, since this is our primary concern in 
the present paper. 



3. Stokes - an overview 

The computer program Stokes performs simulations of radiative 
transfer, including the treatment of polarization, for AGNs and 
related objects. The code is based on the Monte Carlo method 
and follows single photons from their creation inside the source 
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last scattering event 



creation of 
a new photon 




scaCtenno events 



Fig. 1. A photon working its way through the model space. 



region through various scattering processes until they become 
absorbed or manage to escape from the model region (Fig. [U. 
The polarization properties of the model photons are given by 
their stored Stokes vectors. 

Photons are created inside the source regions, which can be 
realized by different geometries. The continuum radiation is nor- 
mally simply defined by the index a of an Fy oc v " power law. 
The Stokes vectors of the photons are initially set to the values 
of completely unpolarized light. 

Various scattering regions can be arranged around the 
sources. The program offers e.g. toroidal, cylindrical, spherical 
or conical shapes. These regions can be filled with free elec- 
trons or dust consisting of "astronomical silicate" and graphite. 
A photon works its way through the model region and gener- 
ally undergoes several scatterings. The emission directions, path 
lengths between scattering events, and the scattering angles are 
computed by Monte Carlo routines based on classical intensity 
distributions. During each scattering event the Stokes vector is 
changed by multiplication with the corresponding Mueller ma- 
trix. For dust scattering, absorption is important, and a large frac- 
tion of the photons never reaches the virtual observer. The rele- 
vant cross sections and matrix elements for dust scattering and 
absorption are computed on the basis of Mie theory applied to 
size distributions of spherical graphite and silicate grains. 

If a photon escapes from the model region, it is registered 
by a web of virtual detectors arranged in a spherical geometry 
around the source. The flux and polarization information of each 
detector is obtained by adding up the Stokes parameters of all 
detected photons. If the model is completely axially symmetric 
these can be azimuthally integrated and, if there is plane sym- 
metry, the top and bottom halves are combined. The object can 
be analyzed in total flux, in polarized flux, percentage of po- 
larization, and the position angle at each viewing angle. The 
light travel time of each photon is also recorded, so it is pos- 
sible to model time-dependent polarization (Gaskell, Shoji, & 
Goosmann, in preparation). 

3.1. Monte Carlo method, photon initialization, and sampling 
the free path length 

Using the Monte Carlo method it is possible to generate a ran- 
dom event x according to a given probability density distribution 



4 



Goosmann & Gaskell: Modeling AGN Polarization 



p(x). Let p(x) be defined on the interval [0, JCmax]- We can then 
construct the probability distribution function P(x) and relate it 
to a random number, r, between and 1 as follows: 



r = P(x)^^ [ p(x')dx'. 
L Jo 



(1) 



The constant C is a normalization constant resulting from 
integration over the whole definition interval [0, Xmax]- Given the 
random number, the corresponding value of x for a single event 
is obtained by inverting equation ([T)- A good description of the 
Monte Carlo method can be found in Cashwell & Everett ( I1959I ). 
In the following, we describe the main routines of Stokes and 
denote all random numbers computed from equation ([T]) by n, 
with / = 1,2,3.... 

To generate a model photon, its initial parameters of po- 
sition, direction of flight, and wavelength all have to be set. 
Different geometries for the continuum region, broad-line re- 
gion, and narrow-line region are available in Stokes. Assuming 
a constant density of the emitting material, a random position 
for the new photon is sampled. The flight direction is given by 
two angles, 6 and 0, defined with respect to a standard polar 
coordinate system. Assuming isotropic emission, the sampling 
equations for the angles are as follows: 



<f> 



arccos(l - 2ri), 
27rr2. 



(2) 
(3) 



The wavelength of the photon is sampled according to the 
intensity spectrum over a range [/Imin, 'imaxl- This leads to: 



A = 



\a". + (a" - A". ]]' , for a 5t 1, 

L mm \ max min/J ' ' 

^minfj^y' , for a = 1. 



(4) 



Here, a denotes the usual power law index of the intensity 
spectrum. 

If we ignore scatterings back into the beam, the intensity of 
a photon beam traversing a slab of scattering material with par- 
ticle number density and cross-section cr will drop by a factor 
of with / being the distance traveled inside the scattering 
region. From this, one can derive the sampling function of /: 



/= -^ln(l -r4). 

Ncr 



(5) 



The factor ^ is the mean free path length. Depending on 
the scattering material, the program uses either a dust extinction 
cross-section o-gxt computed from Mie theory or, in case of elec- 
tron scattering, the Thomson cross section ctes- 

3.2. Polarization formalism, scattering, and photon detection 

The polarization properties of the photons sampled in Stokes 
and their transformation during scattering events rely on previ- 
ous work described e.g. in Fischer, Henning, & Yorke (1994). 
The theoretical basis for the formalism presented in this section 
can be found in Bohren & Huffinan ( I1983I I. 

If we consider a photon being scattered off" a spherical par- 
ticle (see Fig.|2]i, the outgoing electromagnetic wave associated 
with the photon can be resolved into two components, E^^ and 
E_i. These components refer to directions of the electric-field 
vector parallel and perpendicular to the scattering plane. For 



scattering off" a spherical particle, the following relation between 
the incoming and scattered electric fields holds: 



SliO) 
SM 



Eli 



(6) 



The scattering matrix elements, S ](0) and 82(0), are inde- 
pendent of the azimuthal angle 0. In case of Thomson scattering, 
there absolute values obey to simple analytic expressions: 



(7) 
(8) 



For dust scattering, the albedo and the matrix elements of 
a standard dust grain are calculated from Mie theory (see sec- 
tion |33]l. The albedo at the photon wavelength is compared to 
a random number, r^, in order to decide whether the photon is 
absorbed or scattered. If the photon is absorbed it is lost, and the 
cycle starts over with the generation of a new photon. 

The polarization vector of each photon lies perpendicular to 
its trajectory, inside the so-called polarization plane, and denotes 
the preferred direction of the ^-vector It is defined with respect 
to a co-moving coordinate system. The polarization information 
of each photon is coded by four Stokes parameters /, Q, U, and 
V, representing the 4-dimensional Stokes vector We assume that 
newly created photons coming from the source are unpolarized. 
Hence, their Stokes vectors have the simple form: 



(9) 



With each scattering event, the co-moving coordinate system 
undergoes a double rotation: the first rotation, by the azimuthal 
angle <p, occurs around the current flight direction of the pho- 
ton. It rotates the E-vector inside the polarization plane to the 
position of the new scattering plane (see Fig. |2]i. Physically, it 
does not affect the polarization state, although the Stokes vector 
undergoes the following coordinate transformation. 



( ^ 1 




( ^ ] 


Q 







u 







[vj 




loj 



f I* ) 


I 


Q' 




U" 




[v* ) 


V 





COS 2 




sin2(/> 



- sin 2(f> cos 20 







1 









gin 











(10) 



The second rotation occurs in the scattering plane by the 
scattering angle 0. The change of the Stokes vector is determined 
by the Mueller matrix, which for scattering off" a spherical parti- 
cle has the form: 







gout 


1 






[ y™' J 


V 



Su 
Sn 







5 12 

S22 







S 33 



-s. 





S 34 





( ^* 1 




Q* 




u* 




[ y j 



(11) 



The entries of the Mueller matrix are obtained by simple 
relations from the elements of the scattering matrix, S i{0) and 

S2{0)- 

The angle-dependent classical intensity distribution of a 
scattered electromagnetic wave measures the probability of find- 
ing a scattered photon at a given direction. Such probability den- 
sity distributions are derived from equation (|6]l. An important 
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3.3. Computation of dust properties 



z 


y 





—^80° - 




First turn by (j) bringing £^ 
into the new scattering plane 



Fig. 2. Geometry and denotations for a single scattering event. 
The inset shows the first rotation of the ^-vector by the angle <p, 
the view on the polarization plane is along the negative z-axis. 



aspect included in Stokes is that the scattering direction is sam- 
pled depending on the incident polarization vector. The degree 
Pi and the position angle yi of the polarization before scatter- 
ing are computed from the incident Stokes vector and enter the 
sampling equations of and (p - (p' + 180° - yi: 



Piff) = N (|5i(0')l' + |52(0')l')sin(0'W, 



|5i(0)|2 + |52(0)P ' 2 



(12) 
(13) 



The number is a normalization constant in order to have 
( fT2] l range from to 1 for scattering angles between 0° and 180°. 
To sample and 6, the right hand-sides of these equations have 
to be set equal to random numbers r^, r?. The equations are then 
solved for the angles. 

Note that the sampling is independent of the incident polar- 
ization for 6 but not for (p. In several Monte-Carlo polarization 
transfer codes described in the literature, the incident polariza- 
tion does not affect the sampling of the scattering angles. This 
does not present a problem if one considers unpolarized inci- 
dent radiation and low optical depths. Also for very high op- 
tical depths, when multiple-scattering neutralizes the polariza- 
tion inside the scattering region, the incident polarization can be 
neglected. However, results for intermediate optical depths are 
sensitive to the sampling method and they should consider the 
polarization state of the incident photon. 

When a photon escapes from the model region it is recorded 
by one of the virtual detectors. It is then necessary to rotate the 
polarization plane around the flight direction until it matches the 
reference axis of the detector. The Stokes vectors of all incoming 
photons can finally be added up to the values /, Q, U and V. The 
net polarization properties are derived from: 



P = 



r = 



1 u 

— arctan — . 

2 Q 



(14) 
(15) 



Mathis, Rumpl, & Nordsieck ( 119771 MRN) suggested dust com- 
positions to reproduce extinction curves observed in our Galaxy. 
They assumed various types of dust grains having a size distri- 
bution proportional to a', with a being the grain radius and s 
an arbitrary power-law index. Our parameterization of the dust 
properties follows that of MRN and gives a good description of 
observed Galactic extinction curves. The user can choose the ar- 
bitrary minimum and maximum radii of the grain-size distribu- 
tion, its power law index, and the relative abundances of graphite 
and "astronomical silicate". 

The results from Mie scattering theory, i.e., scattering and 
extinction cross sections, albedos, and elements of the scatter- 
ing matrix, are computed using the code given by Bohren & 
Huffman ( 119831 ). We imported complex dielectric functions for 
graphite and silicate measured by Draine & Lee (1984). For 
graphite, two dielectric functions have to be considered since the 
optical properties for light polarized parallel and perpendicular 
to the crystals axis differ from each other. The code therefore 
works with the two different graphite types having abundances 
in a ratio of 1:2. It computes a weighted average for the dust 
composition and grain size distribution defined. The procedure 
is described, for example, in Wolf (1200 3l l. The properties of the 
resulting "standard dust grain" are then used by all dust-related 
routines of Stokes. 

We confine ourselves to using standard Galactic dust such 
as is seen in the solar neighborhood, even though there is evi- 
dence that the tori of AGNs might have different compositions 
and grain size distributions (see Czerny et al. 120041 Gaskell et al. 
2004, and Gaskell & Benker [2006t . Following Wolf & Henning 
j 19991 we parameterize Galactic dust by a mixture of 62.5% car- 
bonaceous dust grains and 37.5% "astronomical silicate". We 
consider grain radii, a, from 0.005 yum to 0.250 fim with a dis- 
tribution n(fl) oc a" with s = -3.5. The resulting cross-sections 
and the albedo are shown in Fig. [3] as a function of wavelength. 
The figure shows that for this particular dust model, the albedo 
is rather flat with a value of 0.55-0.6 over the wavelength range 
considered. For wavelengths < 2500 A it falls to 0.4. The cross- 
sections all decrease regularly with wavelength, with the excep- 
tion of the well-known hump around 2175 A. 

4. Simulation of torus geometries 

In this section we investigate how much of the polarization prop- 
erties of type-1 and type-2 AGNs can be produced by a uniform- 
density torus alone. Kartje ( 119951 ) modeled the polarization in- 
duced by scattering off a cylindrically shaped torus. Their torus 
model was adopted from a fit to NGC 1068 given by Pier & 
Krolik (1992). This torus is geometrically rather compact and 
is located within a radius of 1 pc from the central source. Such 
a cylindrical torus is not necessarily physical, so we examine 
whether the results of Kartje can be confirmed with more gen- 
eral tori, and we extend the range of parameter space explored. 

4. 1 . Curved surfaces versus sliarp edges 

The dusty tori examined by Kartje dl 9951 1, Wolf & Henning 
(fT999b . Young ([2000b . and Watanabe et al. ([2003T l have rather 
sharp edges, and, since we find that polarization results can de- 
pend strongly on geometrical details, we have investigated a less 
artificial torus geometry with an elliptical cross-section. To ex- 
amine the influence of sharp edges of the cylindrical torus on the 
polarization, we define a torus with similar dimensions, and the 
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Fig. 3. Characteristic properties of the dust composition adopted 
for this paper as a function of wavelength. Top: albedo value. 
Bottom: cross-sections for extinction (black, solid), scattering 
(red, dashed), and absorption (blue, dotted). 




Fig. 4. Geometry of the three torus models we consider: (1) 
the cylindrical torus used by Kartje, (2) a compact elliptically- 
shaped torus, and (3) an extended elliptical torus. All tori have 
the same half-opening angle ©o 



same optical depth in the V band (ry) ~ 750 along the radius 
in the equatorial plane. Thus, practically no photon is able to 
penetrate through the torus and only scattering off its surface is 
relevant. However, our torus has an elliptical cross-section (see 
Fig-0 instead of the rectangular cross-section used by Kartje. 

Our results compare very well to those obtained by Kartje 
( I1995I I. In Fig. |5] we show polarization and flux (normalized 
to the flux of the central source) versus wavelength at different 
viewing directions. The torus considered has a half-opening an- 
gle of 00 - 30°. The positive values of P denote that the polar- 
ization vector is oriented perpendicularly to the symmetry axis 
(type-2 polarization). In our simulations the torus is filled with 
standard Galactic dust, parameterized as described at the end of 
section[3] We sample a total of 10^ photons and record spectra at 
10 different viewing angles scaled in cos /, where / is measured 
from the axis of the torus. We show our results as a function of 
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Fig. 5. Modeling a cylindrical torus with an elliptical cross- 
section and 00 - 30° (see section 14. It . Top: polarization, P. 
Bottom: the fraction, F/Ft, of the central flux, F,, seen at dif- 
ferent viewing inclinations, /. Legend: / = 87° (edge-on) (black 
crosses), / = 76° (orange triangles with points down), / - 70° 
(intermediate) (maroon stars), / - 57° (purple triangles with 
points up), i = 41° (green diamonds), / - 32° (red squares), 
and i = 18° (face-on) (blue circles). 



cos / because it gives equal flux per bin for an isotropic source 
located at the center of the model space if there is no scatter- 
ing. Our figure is quite similar to the corresponding diagrams in 
Kartje's paper (see his Fig. 5). 

The only difference between our results and those of Kartje 
is that we generally obtain slightly lower polarization degrees 
and a slightly diff'erent wavelength-dependent slope for the scat- 
tered flux. This can be explained by the fact that we calculate 
our cross-sections from Mie theory of a specific dust composi- 
tion whilst Kartje used cross-sections given by Mezger, Mathis, 
& Panagia ( |T98^ . 

We also investigated the polarization of a compact torus with 
an elliptical cross-section for changing 6o. Again we obtained 
similar results (not shown) to those for Kartje's cylindrically- 
shaped tori. Thus, the differences in polarization between the 
elliptical and cylindrical tori are negligible. Having sharp edges 
in the cylindrical model rather than the more realistic rounded 
edges of the elliptical torus does not introduce spurious effects. 

4.2. The effect of the shape of the inner edge of the torus 

A real torus is undoubtedly thicker than the geometrically thin 
cylindrical torus of Kartje. Direct imaging of NGC 4261 (= 
3C 270) shows that the dusty torus in that AGN extends out 
to 230 pc (Ferrarese, Ford, & Jaffe 1 19961 1. A similar dust lane 
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across the nucleus of M 51 (= NGC 5194) extends by ~ 100 pc 
(Ford et al. '1992'). The inner radii of tori are obtained by infra- 
red reverberation mapping of the hot dust and are in the range 
of tens to hundreds of light-days for Seyfert galaxies (see Glass 
12004 and Suganuma et al.l2006l). 

The outer regions of tori have considerable optical depth, so 
their precise shape is unimportant, since no photons escape par- 
allel to the equatorial plane of the torus. The shape of the inner 
region facing the central energy source is more relevant. Current 
torus models commonly consider inner surfaces that are con- 
vex towards the central source. We thus model optic ally- thick, 
uniform-density tori with elliptical cross sections, an inner ra- 
dius of 0.25 pc, and an outer radius of 100 pc. We compare 
these results to the modeling of a more compact torus with the 
same half-opening angle, Oq - 30°, as in section 14.11 We de- 
termine the dust density by fixing ry at ~ 750. Variability ob- 
servations imply that the size of the optical and UV-continuum 
source in Seyfert galaxies is less than a few light-days, as is also 
expected from simple black-body emissivity arguments. Hence, 
when considering scattering off the torus, we can neglect the fi- 
nite size of the continuum source in our model and assume a 
point-like emission region. Note that this consideration remains 
valid for objects with higher luminosities because both the size 
of the central emission region and the inner radius of the torus 
scale with luminosity. 

The resulting spectra at different incUnations are shown in 
Fig- IS] If the viewing angle, /, is less than Oq (thus corresponding 
to a type-1 object), we only observe a regular type-1 spectrum. 
We find that there is no significant polarization in this case. If we 
look at a type-2 object at a higher inclination angle, only scat- 
tered (and hence polarized) light is detected. This is analogous 
to the results obtained for the compact torus shown in Fig. |5] 
The overall shape of the polarization spectrum for both sizes of 
the torus is rather similar as well. With increasing viewing angle 
the level of the polarization spectrum rises, reaches a maximum, 
and decreases again towards edge-on lines of sight. The shape of 
the P-spectrum does not change significantly between different 
type-2 inclinations. 

There are differences between our results of modeling a large 
torus (case 3 in Fig. HJl with half-opening angle Oq = 30°, and 
the analogous compact torus (case 2 in Fig. |4| with identical 
half-opening angle but smaller dimensions. A striking difference 
occurs in the angular flux distribution: the large torus scatters 
considerably fewer photons towards an observer at intermedi- 
ate viewing angles because they hit the outer parts of it (see the 
illustration in Fig. |2l). Towards edge-on viewing directions the 
probability of seeing scattered photons is much lower than for 
the small torus. The spectral slope of the scattered radiation also 
differs between the two tori. While the spectrum is flat in the case 
of a compact torus it rises towards the blue for the large torus. 
This can be explained by the increasing tendency of forward- 
scattering at shorter wavelengths. Photons escaping at higher in- 
clinations have to undergo back-scattering; this is more likely to 
happen at longer wavelengths. 

There are also differences in the polarization signatures of 
both tori. Although the overall spectral dependence of P is the 
same, the level of P is changed. The strongest changes are at 
higher inclinations when the central source is becoming ob- 
scured by the torus. As with the total flux (see above), for the 
larger torus, P is significantly lower (compare the case of / - 70° 
between the upper panels off Fig. |5]and Fig. |6]l. For a large torus, 
our current models sampling several 10'' photons do not con- 
strain the polarization well at very high inclinations. The number 
of photons scattered into these directions is too small to allow for 
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Fig. 6. Modeling a large torus with an elliptical cross-section and 

00 = 30° (see section 14. 2t . Top: polarization, P. Bottom: the 
fraction, F/Ft, of the central flux, F,, seen at different viewing 
inclinations, Legend: / = 70° (intermediate) (maroon stars), 
/ - 63° (pink triangles with points to the right), / - 57° (purple 
triangles with points up), / = 49° (brown triangles with points to 
the left), / = 41° (green diamonds), / = 32° (red squares), and 

1 - 18° (face-on) (blue circles) 




Fig. 7. Comparison of the compact and the extended torus with 
6*0 = 30° in the V-band. 



sufficient statistics. On the other hand, it clearly follows from our 
computations that the spectral flux at angles / > 76° is reduced 
by a factor of almost ~ 2x10^ with respect to the flux of the 
source. Therefore, the polarized flux at these angles is very low. 

We show the differences in V-band total flux and polarization 
between a large and a small torus in Fig. [8] The top panel shows 
the polarization as a function of the viewing angle, and the bot- 
tom panel the fraction of the light reaching the observer. As was 
shown above, the differences between the two torus shapes are 
most important at higher inclinations. At / ~ 70° the degree of 
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Fig. 8. Differences between large and small tori with an elliptical 
cross-section in the V-band (see section l4~2b . Top: polarization, 
P, and bottom, the fraction, FjFt, of the central flux. Ft, as a 
function of viewing inclinations, /. The dashed lines denote the 
thin elliptical torus (case 2), the solid line the extended torus 
(case 3) 



polarization reaches a difl'erence of 6%, and the flux difi'ers by a 
factor of almost 100. 

4.3. The effect of the torus half-opening angle 

Kartje (1995) has shown that the half-opening angle, do, of the 
torus is an important parameter for the obscuration and reflection 
properties. While modeling large tori, we examine half-opening 
angles ranging from 10° to 75°. Variation of 6o is realized by 
changing the vertical half-axis of the elliptical torus cross sec- 
tion. The other model parameters are defined as for the previous 
case of 6*0 = 30° in section I4r2] 

4.3.1 . Tori with narrow or wide openings 

For large tori, 60 is a dominant parameter for both the degree of 
polarization and the position angle, y. In Fig.|9]we show the po- 
larization of the scattered radiation as a function of wavelength 
and for various Oq. Due to a similar overall shape of the wave- 
length dependence of P, we average the polarization over type-2 
viewing angles, /, with / > Oq. We thereby exclude the high- 
est inclinations with an insuflicient number of photons, where 
the statistics of P are too poor For viewing angles with ; < 60 
(corresponding to type-1 objects seen face-on) the polarization 
is negligible. 
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Fig. 9. Polarization averaged over type-2 viewing angles (see 
section |43] l. A positive value of polarization denotes an /?- vector 
oriented perpendicular to the torus symmetry axis; for negative 
values the /J-vector is aligned with the projected axis. Legend: 
00 = 10° (black dashed Une), = 20° (solid red fine), = 30° 
(green dot-dashed fine), 0o = 45° (blue dots), 6*0 - 50° (long 
yellow dashes), = 60° (brown double dots and dashes), and 
6o = 75° (pink double-dashes and dots). 



Varying the opening angle shows several important things. 
For 6o < 53° the absolute value of the polarization decreases as 
the opening angle increases (see Fig. |9]l, as was found by Kartje 
(11995 ) for compact tori. The polarization vector is oriented per- 
pendicularly to the axis for all viewing directions / > 60, as is 
observed in type-2 AGN. For 60 > 60°, only parallel polariza- 
tion vectors can be seen at viewing angles / > Oq. In this range of 
the absolute degree of polarization increases with the opening 
angle. 

The reason for the flip of the relative position angle can be 
explained by the scattering phase function, and by the geome- 
try of the inner parts of the torus (Kartje 1995). For a distant 
observer looking at the torus along an off'-axis line of sight, the 
scattered radiation comes from the inner surface walls. In part, 
these consist of the inner torus wall facing the observer most 
directly, but they also consist of the two surfaces on the side. 
Due to the scattering geometry, the photons scattered off" the side 
walls are polarized along the projected symmetry axis, whilst the 
photons coming from the far wall are perpendicularly polarized. 
The ratio of the solid angle that the far side of the visible in- 
ner surface subtends to the solid angle that the visible inner side 
walls subtend changes with the half-opening angle of the torus, 
and so does the overall polarization vector. 

4.3.2. Transition case: intermediate torus lialf-opening 
angles 

For intermediate opening angles with 53° < do < 60° the orien- 
tation of the polarization position angle seen at type-2 viewing 
angles depends on the exact inclination. We illustrate such a case 
in Fig. [To] where we set 60 = 57°. 
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Fig. 10. Modeling an expanded torus at an intermediate opening 
angle of Oq = 57°. The graph shows the polarization degree in 
the visual band versus inclination angle /. A positive value of 
polarization denotes an E-vector oriented perpendicular to the 
torus symmetry axis; for negative values the E-vector is aligned 
with the projected axis. 



For a line of sight passing close enough to the horizon of 
the torus (i.e., when / is only moderately larger than ^o) we find 
that the polarization vector is parallel, which means that type- 
1 polarization can be produced at obscured viewing inclinations 
(Fig.fTOb. If the inclination increases further the polarization vec- 
tor switches back to type-2 polarization. It is interesting to note 
that such a torus can produce significant polarization degrees up 
to 2% for both orientations of the £- vector. 

In order to illustrate the integral effect of the opening angle 
on the polarization, we plot in Fig.[TT]the polarization, feft, av- 
eraged over all type-2 viewing positions, and over wavelength, 
as a function of the half-opening angle of the torus. The differ- 
ence between type-1 and type-2 polarization is ignored in Peff. 
The absolute values of P are integrated. 

The figure shows that the torus polarizes most effectively 
when having either a small or a large half-opening angle. In the 
transition region between type-1 and type-2 polarization (i.e. for 
53° < 9 < 60°) the integrated polarization goes through a mini- 
mum. 

4.4. Wavelength insensitivity of polarization due to dust 
scattering 

Wavelength-independent polarization is widely taken to be a sig- 
nature of electron scattering, but we have shown in Figs. 6 and 
9 that dust scattering can also produce wavelength-independent 
scattering. Thus a flat polarization curve is not a unique signa- 
ture of electron scattering. As Zubko & Laor (2000) point out, 
the wavelength dependence of polarization provides a probe of 
the grain scattering properties. Inspection of the wavelength- 
dependent polarization curves for the large torus geometries con- 
sidered above (see Figs. 6 and 9) shows that the polarization 
for half-opening angles with 30° < Oq < 60° is wavelength- 
independent over the optical and most of the UV. For values of 
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Fig. 11. Effective polarization, Pgif (see section [43] ), for type-2 
viewing positions as a function of the half-opening angle of an 
extended torus (case 3). 



Oo outside this interval the wavelength-dependence of P is rather 
low and does not exceed a factor of 2. Since the scattering cross 
section of interstellar grains increases strongly from the optical 
to the UV, wavelength-independent polarization is commonly 
supposed to be the fingerprint of electron scattering. However, 
scattering in opaque dust clouds produces relatively grey scat- 
tering (Kishimoto 2001). 

Our apparently contradictory result of relatively wavelength 
independent polarization with dust scattering arises because we 
are considering scattering off optically-thick material, and be- 
cause of the relatively small variation of the albedo over the 
optical and UV spectral regions (see Fig. [3]). The approximate 
constancy of the albedo is because the scattering and absorption 
cross-sections vary in a similar manner with wavelength. Since 
we assume an optically-thick torus, we see emergent photons 
that have been scattered at an optical depth r ~ 1 . This is regard- 
less of wavelengtlQ. The increase in scattering cross-section with 
decreasing wavelength only means that the shorter wavelength 
photons we see have been scattered closer to the surface of the 
torus. 

A significant change in albedo with wavelength, however, 
will cause a color dependency in the intensity and polarization 
of the scattered ligh£|. Shortwards of ~ 2500 A the albedo de- 
creases, but this range is at the lower limit of the spectral range 
considered in our modeling. The effect can be seen in the nor- 
malized flux spectra of the torus models shown in Figs. |5] and 
|6] For the Galactic dust composition we implemented, it is less 
visible in the polarization spectra. 

Another grain property that needs to be considered is the 
degree of asymmetry of the scattering since this is effectively 
an angle-dependent albedo change. Toward shorter wavelengths. 
Galactic dust grains are more strongly forward scattering and 



- This is the reason that the sunht sides of clouds in the earth's atmo- 
sphere are extremely white. 

^ This is the cause of colorations in the atmosphere of the giant plan- 
ets. 
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the polarization phase function changes (see Draine |2003l l. As 
for Thomson scattering, forward-scattered light has a lower po- 
larization than sideways-scattered light. The polarization spectra 
obtained for the torus models depend on these phase functions. 
They additionally explain why a slight wavelength-dependence 
of the polarization can be found for very narrow or very wide 
opening angles of a large torus (see Fig.|9l). 

5. Polarization from polar-scattering regions 

Scattering in polar regions of AGNs has allowed the discovery 
of hidden Seyfert-1 nuclei in type-2 objects by radiation being 
periscopically scattered around the obscuring torus. The central 
parts of the polar double cone have to be at least moderately 
ionized due to the intense radiation from the AGN. The medium 
could be associated with the warm absorber seen in many AGN 
(see Komossa 1999 for a review). The Doppler shift of the X- 
ray absorption lines indicates that the medium is outflowing at 
roughly 1000 km/s. With increasing distance from the center, the 
outflow velocity and intensity of the radiation decrease. Beyond 
the sublimation radius, dust could also be present. However, this 
dust must be optically thin, as type-1 objects are not obscured. 

5. 1. Polar electron scattering 

The polarization induced by scattering in polar, conical electron- 
scattering regions has been the subject of several previous stud- 
ies. Brown & McLean (! 19771 1 developed a formalism to compute 
the polarization expected from scattering inside optically thin, 
axisymmetric scattering regions. This formalism was applied by 
Miller & Goodrich ( 1 990) and MiUer et al. ( 1 99 1 ) to compute the 
polarization for polar scattering cones. Wolf & Henning (1999) 
and Watanabe at al. (2003) extended the modeling to optically- 
thick material using Monte-Carlo techniques that can account 
for multiple scatterings. 

We confirm such results in Fig.[T2]using Stokes. The figure 
shows the degree of polarization and the total flux as a function 
of the observer's inclination for an electron scattering double- 
cone of uniform density and with the optical depth Tes = 1 . This 
optical depth is measured in the vertical direction between the 
inner and the outer shell of one cone. In order to isolate the ef- 
fects of the scattering cone from the polarization induced by the 
disk, we use an anisotropically emitting central source with the 
emission angles being restricted to the solid angle defined by 
the scattering cones. The three curves denote the half-opening 
angles 0c = 10°, 0c = 30°, and Oc = 45°. The inclination is 
measured from the symmetry axis of the double-cone. 

As expected, polar electron-scattering cones produce type- 
2 polarization directed perpendicularly to their symmetry axis. 
The degree of polarization rises from face-on to edge-on view- 
ing angles. The latter effect is due to the angle-dependent polar- 
ization phase function of Thomson scattering. For wider open- 
ing angles of the cones, the net polarization P decreases because 
it is the result of integrating a broader distribution of polariza- 
tion vectors. The breaks of the polarization curves at ; = 6c 
in Fig. [12] are due to the impact of multi-scattering inside the 
cones, the analogous breaks in total flux curves are due to the 
angle-restricted central emission. 

In Fig.[T3]we plot the influence of the optical depth on the po- 
larization for the polar-cones with 0c = 30°. The various curves 
denote different optical depths. A similar case was considered 
by Watanabe et al. (120031) . The density of their electron cones 
varies with the distance from the center according to a power 
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Fig. 12. Modeling double polar cones of various half-opening 
angles. Top: polarization, P, with positive values denoting type- 
2 polarization (perpendicular to the symmetry axis). Bottom: the 
fraction, F/Ft, of the central flux. Both values are plotted versus 
the inclination / with respect to the observer The different sym- 
bols denote diff'erent half-opening angles of the double-cone. 
Legend: 6c = 10° (black circles), 6c - 30° (red squares), and 
6c - 45° (blue diamonds). The optical depth between the inner 
and outer shell of the cones is set to Tes = 1 . 



law. Comparison of Fig. [T3] with Fig. 7 of Watanabe et al. (2003] 
bottom panel) shows that the diff'erence in P is very small with 
respect to a uniform-density torus. 

5.2. Polar dust scattering 

Beyond the dust sublimation radius the scattering cone could 
contain dust. We investigated this using a similar bi-conical ge- 
ometry and our Galactic dust prescription - see section [3] In 
Fig. [14] we show the polarization and total flux resulting spec- 
tra for a centrally-illuminated dust cone seen at different view- 
ing angles /. The half-opening angle of the cone has been set 
to 6c = 30°, and its optical depth in the V-band to the moderate 
value of Tv = 0.3. The strong wavelength dependence of the dust 
extinction properties (see Fig.[3]l is clearly visible in the figure. 
It differs for polar viewing angles, which cross the cone, from 
those along equatorial directions. The former ones show the dust 
extinction seen in transmission, while the latter ones show dust 
reflection. 

The total flux shown in the bottom panel of Fig.[T4]is signifi- 
cantly reddened when i < 6c. In addition to that, the well known 
extinction feature at 2175/im is seen. Its depth decreases with in- 
creasing inclination. In reflection (i.e., along equatorial viewing 
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Fig. 13. Polarization degree by polar electron scattering cones 
with half-opening angle 0c - 30° plotted versus the inclination 
; with respect to the observer The positive values denote type-2 
polarization (perpendicular to the symmetry axis). The different 
symbols denote different optical depths between the inner and 
the outer radius of the cone. Legend: tes = 0.01 (black circles), 
tes - 0.3 (red squares), tes = 1 (blue diamonds), tes - 3 
(green triangles with points up), and tes = 10 (triangles with 
points left). 



angles) the spectra are quite different. They peak at ~ 2500 A, 
where the albedo goes through a maximum (see Fig. [3]), and 
level down slowly toward longer wavelengths. This behavior is 
the same at all inclinations with / > 9c\ the spectra differ just 
slightly in normalization. 

The polarization spectra remain below 1 % when the double- 
cone is seen in transmission because forward scattering does not 
induce significant polarization. In reflection, however, polariza- 
tion becomes significant, and, for all inclinations with / > 0c, 
it rises towards longer wavelengths. The shapes of the spectra 
are, again, not very dependent on the inclination; they only dif- 
fer in normalization. The highest level of polarization is obtained 
when the cone is seen edge-on, which corresponds to perpendic- 
ular scattering angles. 

Note that, while the wavelength dependence in polarization 
and flux is different from the one obtained for the dusty torus in- 
vestigated in section|4](because we are considering optically thin 
material in the cones), the of the flux rise to shorter wavelengths 
is still less than the rise in the cross sections. 



6. Polarization from equatorial electron 
distributions 

The major difficulty in modeling the type-l/type-2 polarization 
dichotomy is producing polarization parallel to the symmetry 
axis. We have illustrated above that polar scattering cones can 
only produce polarization vectors oriented perpendicular to the 
axis, and that the polarization produced by the dusty torus at 
type-1 viewing angles is mostly negligible. An exception was 
the case of extremely-thin tori, but these are inconsistent with 
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Fig. 14. Modeling a dusty double cone of half-opening angle 
0c - 30°. Top: polarization, P. Bottom: the fraction, FjFt, of the 
central flux, F,, seen at different viewing inclinations, /. Legend: 
/ - 87° (edge-on, black crosses), / - 76° (orange pluses), / = 63° 
(pink triangles with points to the right), / - 49° (brown triangles 
with points to left left), / = 32° (red squares), and i - 18° (face- 
on, blue circles). 



the observations for at least two reasons: they are incompatible 
with the ratios of type-1 and type-2 objects, and very flat distri- 
butions have too low covering factors to produce the reprocessed 
IR emission. 

It is therefore necessary to introduce a third type of scattering 
region. It was pointed out a long time ago that parallel polariza- 
tion can be produced by a thin emitting and scattering disk (e.g., 
Chandras ekhar . 1 960. Angel , 1969 Antonucci 1983 , Sunyaev & 
Titarchuk ri9851 l. Goodrich & Miller ( |1994| i suggested that the 
parallel polarization of type-1 objects arises from scattering in a 
flattened equatorial medium located around the accretion disk or 
even the BLR. This model has been investigated in a series of pa- 
pers by Young .2000.. Smith et al. 2 002. .2004. and l2005T l using the 
GSM. In this model a rotating equatorial scattering disk explains 
the velocity-dependent polarization structure across the broad 
emission lines of a large sample of Seyfert galaxies very well, 
and also helps explaining the type-l/type-2 dichotomy. Unlike 
the GSM, our code includes the effects of multiple scattering, so 
we can consider scattering regions with higher optical depths. 
We first constrain possible contributions of an AGN accretion 
disk to the polarization (16. It . Then, we investigate with Stokes 
the effect of equatorial scattering disks on the continuum radia- 
tion (section |672] l. 
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Fig. 15. Polarization versus viewing angle for a geometrically- 
thick emitting and scattering electron disk with ^ = 0.5. The dif- 
ferent curves denote various vertical optical depths for the disk. 
Positive polarization values stand for polarization position an- 
gles perpendicular to the disk's symmetry axis, for negative val- 
ues the polarization vector is aligned with this direction. Legend: 
tes - 0.05 (black circles), tes - 0.5 (red triangles with points 
up), Tes = 1-25 (blue squares), tes - 1 (green triangles with 
points left), tes = 1.5 (orange diamonds), tes = 2.5 (brown tri- 
angles with points down), tes - 5 (purple crosses), tes = 10 
(maroon stars), and tes = 25 (cyan triangles with points right). 



6.1. Emitting and scattering disks 

We consider emitting and scattering disks with a ratio ^ of to- 
tal disk height h to diameter d. The disk is modeled in a plane- 
parallel approximation. It has a cylindrical shape and is unifor- 
mally filled with electrons. We consider a "thick disk" cross sec- 
tion with J = 0.5 (see Fig.[T5]l and a "thin disk" with ^ = 0.01 
(see Fig. [16T l. The various curves in the diagrams refer to dif- 
ferent vertical electron-scattering optical depths, tes, measured 
vertically from the central plane to the surface. The optical depth 
is varied by adjusting the electron density. We investigate Tes be- 
tween 0.001 and 50. 

For the case of a geometrically-thick emitting disk (see Fig. 
[TSl ) we find moderate polarization values of at most a few per- 
cent. The degree of polarization depends strongly on the view- 
ing direction and the optical depth. For lower optical depth the 
vector is aligned with the disk's symmetry axis. It flips to a 
perpendicular orientation for tes greater ~ 10. 

Much stronger values of the polarization can be obtained 
when the disk is flatter For the geometrically-thin disk we have 
a similar qualitative behavior as for the thick disk (compare 
Fig. [16] with Fig. [T5] l, but the degree of polarization reaches 
higher values and is significant even for near to face-on viewing 
directions. The flip to perpendicular (type-2) polarization occurs 
at a moderate optical depth (tes ~ 5) and at an edge-on viewing 
angle that cannot be seen for type-1 objects. A thin disk with 
moderate optical depth will thus produce parallel polarization 
for all type-1 viewing positions. 
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Fig. 16. Same parameters as in Fig. [15] but for a thin disk with 
^ = 0.01. Legend: tes = 0.001 (black circles), tes = 0.01 (red 
triangles with points up), tes = 0.05 (blue squares), tes = 0.2 
(green triangles with points left), tes = 0.5 (orange diamonds), 
Tes = 1 (brown triangles with points down), tes = 3 (purple 
crosses), tes = 5 (maroon stars), and tes = 10 (cyan triangles 
with points right). 



The polarization behavior of our uniformly emitting disks 
can be explained in the same manner as the behavior of the po- 
larization of the oblate spheroids examined by Angel (1 19691 1. For 
low Tes the net polarization is mainly determined by the photons 
traveling parallel to the disk plane and then being scattered to- 
wards the surface. For an observer, who does not observe the 
disk exactly face-on, the integrated scattered flux from the disk 
surface will be polarized along the projected direction of the disk 
axis. This can be understood by the fact that polarization by elec- 
tron scattering is most eflicient for orthogonal scattering angles. 
For the same reason, the polarization is also strongest at edge-on 
viewing angles. 

When the optical depth becomes higher, multiple scattering 
of photons traveling toward the disk surface becomes relevant. 
The polarization vector induced by the last scattering event be- 
fore leaving the disk will preferably be oriented perpendicular to 
the disk axis. Hence, the polarization position angle flips and on 
the way to this transition P becomes very low. 

The emission and scattering disks investigated in this sec- 
tion are unlikely the cause of type-1 polarization in AGN. 
High enough polarization parallel to the disk symmetry axis is 
only produced in a disk which is geometrically and optically 
thin. While accretion disks according to the standard model of 
Shakura & Sunyaev (1973) are indeed geometrically thin, they 
are, however, optically thick. They produce too little polarization 
at type-1 viewing angles, as one can tell from the curves in the 
left part of Fig. [16] At type-2 angles the polarization can reach 
up to the well-known limit of ~ 11.7% for the highest optical 
depths (Chandrasekh a n 19601 1, but with a perpendicular orienta- 
tion of the polarization vector 
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Fig. 17. The equatorial scattering geometries investigated: (a) 
flared disks, (b) tori, and (c) equatorial cylinders. The definition 
of the half-opening angle Or, is shown for each case. 



6.2. Equatorial scattering wedges 

Instead of scattering from the accretion disk, it has been pro- 
posed that type-1 polarization is caused by scattering by an equa- 
torial wedge (Goodrich & Miller 1994 ) or by a flared equatorial 
disk (Smith et al. .200 2). Here, we use Stokes to model such 
equatorial regions for different geometrical shapes and optical 
depths. 

6.2.1. Flared disks 

The geometrical shape of a flared disk of half-opening angle 6d 
is given in Fig.[T7](a). We assume that the scattering is only off 
electrons. This is a good assumption for AGNs if the scattering 
region is located close enough to the center for hydrogen to be 
mostly ionized. We fix the Thomson optical depth in the equa- 
torial plane between the inner and the outer shell, tes, to have 
tes = 1- 

In Fig. [18] we plot the total flux and polarization versus in- 
clination, for flared disks with various half-opening angles. 
Both the total flux and the polarization depend significantly on 
the half-opening angle of the flared disk, but the polarization is 
always parallel to the symmetry axis of the system. The flared 
disk hence produces type-1 polarization, while the polarization 
degree obtained at type- 1 viewing angles remains rather moder- 
ate at around 2%. 

For viewing angles / < 90° -6u, the central source is directly 
visible, which leads to higher fluxes than towards edge-on view- 
ing angles / > 90° -0d- If the scattering wedge crosses the line of 
sight, a fraction of the radiation is scattered out of the way. The 
polarization degree increases as the inclination goes from face- 
on towards edge-on viewing angles. At face-on inclinations, the 
line of sight is nearly aligned with the axis of the wedge and 
therefore the observer sees a more axisymmetric system. 

For a type-1 viewing direction with / ~ 30° the maximum 
polarization is obtained for values of 6d in the range of 20° - 
30°. It turns out that changing the optical depth from r ~ 1 de- 
creases the polarization. To illustrate this, we plot in Fig. [T9] the 
polarization degree for a flared disk with 6^ -25° versus the 
inclination / for various optical depths. 

As in Fig. [TS] (top), the polarization curve shows a disconti- 
nuity when the line-of-sight passes the horizon of the flared disk. 
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Fig. 18. Modeling equatorial flared scattering disks of various 
half-opening angles 0£, (see section l6.2.1l ). Top: polarization, P, 
with negative values denoting type-1 polarization (parallel to the 
symmetry axis). Bottom: the fraction, F/F,, of the central flux. 
Ft, plotted versus inclination, /, with respect to the observer 
The diff'erent symbols denote different half-opening angles of 
the flared disk. Legend: Or, - 5° (black circles), 9d = 10° (red 
squares). Or, - 20° (green diamonds). Or, - 30° (blue triangles 
with points up), Oo - 40° (orange triangles with points to the 
left). Or, - 50° (purple triangles with points down), and Od - 60° 
(brown triangles with points to the right). 

However, with larger optical depths, the jump in polarization be- 
comes shallower because multiple scattering within the disk has 
a depolarizing effect. For very low optical depths, r < 0.1, the 
scattering disk is practically invisible, as one would expect. 

The flared disk has somewhat similar polarization properties 
to the emitting and scattering disks investigated in section 16.11 
The irradiation pattern is different but is again axisymmetric. 
Therefore, at low and moderate optical depths, the polarization is 
again mainly defined by photons traveling in a direction parallel 
to the disk plane being scattered at orthogonal scattering angles. 
As for emitting and scattering disks, this produces type-1 polar- 
ization. The central funnel of the flared disk presents a similar 
geometry to the dusty tori discussed in section 143] For flat tori 
this produces again type-1 polarization. Because of the funnel, 
the source is directly visible at type-1 viewing angles, which has 
a strong depolarizing effect. This explains, why the net polariza- 
tion of flared disks at type-1 inclinations remains rather smaU. 

6.2.2. Other geometries 

Finally, we investigate how the geometry of the equatorial 
scattering regions influences the polarization for the shapes in 



14 



Goosmann & Gaskell: Modeling AGN Polarization 



0.0 i-e 



-0.02 - 



-0.06 - 



-0.08 - 




J 1 I I L 



_] I L. 



J I L 



10 20 30 40 50 60 70 80 90 



0.03 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 



0.02 - 



■ / 
/ — 



0.01 



- 



-0.01 - 



-0.02„ 



I I I I I I I I I I I I I I I I I I I I I I I I i_ 

10 20 30 40 50 60 



Inclination 1 [deg] 

Fig. 19. The degree of polarization by electron scattering by an 
equatorial flared disk with half-opening angle 6u = 25° plotted 
versus the inclination, /, with respect to the observer. The neg- 
ative values denote type-1 polarization (parallel to the symme- 
try axis). The diff'erent symbols denote different optical depths 
between the inner and the outer radius of the disk. Legend: 
tes =0.1 (black circles), tes = 0.4 (red squares), tes - 0.7 
(green diamonds), tes = 1.0 (orange triangles with points up), 
Tes = 1 .4 (blue triangles with points to the left), tes = 3 (maroon 
triangles with points to the right), and tes = 5 (brown crosses) 



Fig. [T71 For all three shapes the optical depth is defined hori- 
zontally in the equatorial plane between the inner and the outer 
radius of the scattering region, and fixed at tes = 1 . The inner 
and outer radii are equal for the three diff'erent geometries, as 
are the half-opening angles of 0d = 30°. The inclination of the 
system is assumed to be at the maximal type-1 viewing angle of 
~ 30°. 

Fig. l20l shows the net polarization versus 0u, defined for the 
individual shapes as shown in Fig. [IT] It can be seen that except 
for thick flared disks at large opening angles, there is a similar 
dependence of P on the half-opening angle. Changing the geo- 
metrical shape of the scattering disk shifts the overall level of 
polarization somewhat, but the dependence on the opening an- 
gle largely remains the same. In all cases the scattering regions 
produce a polarization vector aligned with the axis of symmetry 
at type-1 viewing angles. 

We have only considered the possibility of electron scatter- 
ing close to the black hole, and not dust scattering, for several 
reasons. The main reason is that dust will not survive this close 
to the central continuum source. IR reverberation mapping puts 
the inner edge of the dusty torus just outside the broad-line re- 
gion in agreement with the predicted sublimation temperature 
(Suganuma et al. 2006). If we replace the electrons in the models 
discussed in this section with dust it does not give a high enough 
polarization (less than around 1 % even for highly-inclined view- 
ing positions and around 0.1% at best for more likely near face- 
on, type-1 viewing angles) because the dust is forward scatter- 
ing. Since there is absorption for dust grains the number of pho- 
tons coming out is also lower than in the electron-scattering case. 



Half-opening angle [deg] 

Fig. 20. The eff'ect of the half-opening angle of the equato- 
rial electron scattering region on polarization for three differ- 
ent geometries viewed from the maximal type-1 viewing posi- 
tion. Legend: flared disk (blue solid line), scattering cylindrical 
disk (dashed red line), and scattering torus (green dotted line). A 
negative polarization denotes polarization parallel to the axis of 
symmetry. 

In summary, the polarization arising from equatorial scatter- 
ing regions has the correct orientation of the E-vector for type-1 
objects, but the degree of polarization remains moderate at type- 
1 viewing angles. This is largely independent of the geometry 
chosen. The low polarizations arise partly because even at max- 
imal type-1 viewing directions the system is still close to be- 
ing axisymmetric. The resulting low polarization is further de- 
creased by dilution with the unpolarized light coming from the 
central source. Multiple scattering will also play a role in depo- 
larizing the radiation. We obtained maximum polarization values 
for Tes ~ 1- 

7. Discussion 

7.1. Polarization from a dusty, optically-thick torus 

From our modeling of uniform-density tori we obtain polariza- 
tions of a few percent from the torus alone for type-2 objects 
(i.e., without the need to evoke electron scattering cones - al- 
though such cones are clearly present in NGC 1068). The wave- 
length dependence of polarization is rather similar for different 
geometrical shapes of the inner torus surface. A flat surface (case 
1), a slightly curved shape (case 2), or a convex shape (case 3) 
all reveal a nearly flat polarization spectrum over the optical and 
UV range for a broad range of half-opening angles. Thus, po- 
larimetry alone does not give the geometry of dust distributions 
in the central regions of AGN. 

We have found that the polarization is somewhat higher for a 
smaller steeper torus (see Fig.[8]l than for a larger more gradual 
one. At the same time, the obscuration efficiency at type-2 view- 
ing angles is much higher for gradually curved tori. A real torus 
probably has a steep inside to it. This is because the dust radius is 
set by dust sublimation and so will be determined primarily by 
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the inverse square of the distance from the central source. The 
inside of the torus is thus probably concave towards the central 
source. The source creates a spherical region with a radius of the 
order of the sublimation radius of the dust grains, where the dust 
composition is likely to be very different to the dust farther away. 
Destruction of small dust grains will modify the grain size distri- 
bution and hence the albedo and scattering properties are likely 
to change. It will be an interesting future task to model polariza- 
tion by scattering off a dusty torus with a more general geometry 
and using dust compositions different from the Galactic compo- 
sition we have adopted here. 

We have constrained the relation between the polarization 
from the torus and the opening angle (see Fig. |9l). For narrow 
opening-angles, the polarization is stronger and has a polariza- 
tion vector perpendicular to the symmetry axis. With increasing 
half-opening angle the polarization decreases and then gradually 
switches to a parallel polarization vector Recent derivations for 
the number ratio of Seyfert-1 to Seyfert-2 galaxies constrain the 
half-opening angle of the torus to a range of 38° - 48° depending 
on the sample and the analysis considered (see e.g., Tovmassian 
[2001 , Schmitt et al. 2001 ). Thus, for Seyfert galaxies, a realistic 
dusty torus can only produce type-2 polarization. 

The fraction of type-2 AGNs decreases with luminosity, de- 
clining from 80% when log L = 42 to 30% when log L = 46. (see 
Ueda et al. 120031 Hasinger 2004 for recent counts). Therefore, 
Lawrence ( |1991| ) proposed the receding torus model, in which 
higher luminosities of AGNs correspond to larger torus half- 
opening angles Oo. This leads to the prediction that, on average, 
the polarization of Seyfert-2 galaxies decreases with luminosity. 
The situation will be complicated by possible additional polar 
scattering that enhances type-2 polarization, however, there has 
to be a minimum amount of polarization coming from the torus 
alone and we predict that this decreases with luminosity. The 
case is different for higher luminosity quasars because there is 
good evidence that they have larger torus half-opening angles 
(see e.g. Simpson 2003). We have shown that if the torus is large 
and elliptically shaped, type-2 quasars with Oq > 60° may show 
polarization vectors aligned with the symmetry axis. Again, this 
idea only holds if polar scattering is low, which does not seem to 
be the case for recently observed of type-2 quasars (Zakamska 
et al. I2005I I. But, in any case, the torus has to contribute to the 
observed polarization to some extent, and the amount of its po- 
larization must increase with the torus opening angle. A thin, 
disk-like torus was considered for the high-polarization quasar 
OI 287 (Rudy & Schmidt [T9881 ). 

A larger opening angle of the torus corresponds to a smaller 
covering factor of the central source with dust. Therefore, one 
measure of the opening angle is the relative IR flux since it de- 
pends on the dust covering factor We predict that the degree of 
polarization, on average, should be correlated with the relative 
strength of the thermal IR flux. So far, R Smith et al. ( 120031 1 
found a correlation between the K-band luminosity and the 
broad band optical polarization of a sample of 2MASS QSOs. 
Their Fig. 4 shows that QSOs with higher K-band luminosity 
have also higher polarization percentages in the optical and the 
authors point out that this could be due to more luminous or 
more extended scattering regions. Polarization due to scattering 
by a dusty torus can explain this correlation for type-2 objects: 
the higher K-band flux corresponds to a larger covering factor 
of the source and thus to a more narrow opening angle of the 
torus. Such tori produce higher degrees of type-2 polarization. 
In the data shown by R Smith et al. (2003) most of the high po- 
larization quasars are indeed of the spectral type 1.5 or 1.8-2. 
However, there are also a few type-1 objects that have a high 



K-band flux. A larger statistical sample of AGN needs to be ex- 
amined for correlations between the optical polarization and the 
IR flux in order to draw a solid conclusion on this issue. 

The lack of wavelength-dependent polarization we find is 
important for interpreting spectropolarimetry of type-2 AGNs. 
Wavelength-independent polarization is observed in the nuclei 
of the Seyfert-2 galaxy NGC 1068 (Antonucci & Mifler, 1985j, 
for example. The wavelength independence of the polarization 
close to the nucleus of NGC 1068 suggests electron scattering in 
a polar mirror (Miller, Goodrich, & Mathews, '1991'; Kinney et 
al. ( 1991 ) and even enabled a 3D-decomposition of the mirror in 
single electron blobs (Kishimoto 1 1 999l l. However, it is difficult 
to maintain the ionization of the outflow, which requires an in- 
trinsically anisotropic ionizing continuum (Miller, Goodrich, & 
Mathews, 1991). We suggest that wavelength-independent po- 
larization also partly arises from the walls of the torus and com- 
bines with the polar scattering. The fact that dust scattering can 
mimic wavelength-independent polarization has been discussed 
for distant radio galaxies (Dey et al. |1996t Cimatti et al. 11 9961 1. 
while it was stressed that a sharp drop of the dust albedo below 
2500 A should lead to a break in the polarization spectrum. In 
our modeling, such a feature is barely visible with the Galactic 
dust composition we used. However, the grain model we applied 
is rather simple in that it only assumes spherical grains. Most 
dust models of this type predict a sharp drop of the albedo be- 
low 2500 A while observationally constrained dust albedos rise 
toward shorter wavelengths in the Far-UV (Gordon 12004! ) . The 
ambiguity between electron scattering and dust scattering for the 
wavelength-dependence of polarization can hence persist even 
below 2500 A. 

7.2. Polarization from polar scattering cones 

We have modeled polar scattering regions considering of both 
dust and electrons. As expected, only type-2 polarization is ob- 
tained, and the polarization degree rises towards edge-on view- 
ing angles. Larger opening angles of the cone lead to lower po- 
larization values. However, the results obtained here should still 
be used with care for the interpretation of real data. It is an im- 
plicit assumption of our modeling (and of previous modeling) 
that the medium is homogeneous and extends over the whole 
volume of the cone. This does not need to be the case as was 
shown for NGC 1068 in which the ionization cones are only 
party filled with scattering material (Capetti 1995a). For cones 
with larger opening angles, the detailed distribution of the ma- 
terial and the fact that it might be clumpy plays a role in de- 
termining the net polarization and simulations should be more 
detailed (Kishimoto I1996I I. For instance. Miller, Goodrich, & 
Mathews (|19911 l observe strongly wavelength-dependent polar- 
ization from a cloud near the nucleus. The polarization rises to- 
ward shorter wavelengths. This is frequently explained by dust 
scattering because the dust scattering cross section increases to- 
ward the blue as well. However, our computations for scattering 
off optic ally- thin dust shown in Fig. [T4l re veal a different behav- 
ior. At wavelengths longer than ~ 2500 A the polarization degree 
seen at type-2 viewing angles decreases toward the blue because 
of the wavelength dependent polarization phase function (Zubko 
& Laor 2000). The reason that the individual cloud in NGC 1068 
shows rising polarization to shorter wavelengths must therefore 
rely on dilution by an unpolarized blue starlight spectrum from 
the host galaxy. 

Another important reason why the polar scattering regions 
should not be homogeneous is the resulting net polarization at 
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type-2 viewing angles. If the cones have a low optical depth, the 
expected polarization would be very high as we show in Fig.[T3] 
This picture can be accurate only for localized scattering regions 
such as the individual clouds resolved in NGC 1068 which have 
a polarization up to 60% (see Fig. 4 in Capetti .l995aJ . For distant 
objects, where the polar regions cannot be resolved, the observed 
net polarization is lower Therefore, as Kishimoto (1996^ pointed 
out, multiple-scattering in a clumpy medium should reduce the 
net polarization. 

We have found in our models that the radial dependence of 
the scattering material inside the cones is not very important in 
determining the net polarization. The important parameter is the 
optical depth. This is true if one considers isolated cones with no 
other scattering mirrors involved. If one also assumes that there 
is an obscuring and reflecting torus, the situation must change. 
The torus will have a collimating effect, particularly on the cen- 
tral parts of the cones, so their central density becomes more 
important for the net polarization. 

Our study of dust scattering in polar cones again under- 
lines the importance of the optical depth for the wavelength de- 
pendence of the resulting spectra in polarization and total flux. 
While an optically-thick obscuring torus produces polarization 
that is roughly constant over the optical and UV range, scattering 
in cones induces a characteristic wavelength dependence. This 
dependence should vary with the dust composition, and future 
observations and modeling could constrain the dust composition. 

7.3. Relative importance of toroidal and polar scattering 
regions 

To a first approximation, the net polarization, P„e,, due to scatter- 
ing in polar cones and scattering off" the torus is given by the sum 
of the polarized fluxes Pcones x F cones and Ptoms x Ftoms coming 
from the two scattering regions. This summation neglects radia- 
tive coupling between the two scattering regions, but if we take 
the polar cone to be optically thin, it should give a good estimate. 
By comparing the two single polarized fluxes, we can constrain 
the optical depth of the polar scattering region at which the torus 
becomes important for the polarization. We assume a Seyfert-2 
galaxy with a torus half-opening angle of - 30° seen at the 
inclination ; = 63°. The polarized flux for the given parameters 
of a large torus can be derived as shown in section |4] We have 
used the same computational method as in section |5] for the po- 
lar cones. The equatorial obscuration is realized by limiting the 
emission angle to 9^). We find that a polar double-cone with the 
same opening angle, 6q - 30°, as the torus can have an optical 
depth as low as tes ~ 5 X 10"^ to produce a comparable polarized 
flux along / - 63°. 

We also give an estimate for the optical depth of the scatter- 
ing cone needed in a Seyfert-2 system with half-opening angle 
00 = 60° seen at an inclination of / - 70°. Again modeling the 
polarized fluxes of both components separately, we need an op- 
tical depth of Tes ~ 5 x 10"^ for the scattering cone to produce 
a comparable polarization to the torus. 

Our modeling thus shows that it is not possible to indepen- 
dently distinguish between polarization by the torus and by polar 
scattering. As soon as there is even very optically thin polar scat- 
tering material, it dominates the net flux and polarization. The 
very efficient obscuration of large tori stresses the importance of 
polar scattering regions for the detectability of type-2 objects, 
and at the same time it lowers the importance of the exact po- 
larization signature at very high inclinations. The polarized flux 
scattered by a torus at extreme edge-on viewing angles is so low 
that it hardly contributes the net polarization at all. 



Both types of scattering regions are closely related, however, 
because the torus collimates the light entering the polar cones. 
Since the albedo of Galactic dust ranges between 40% and 60% 
a significant fraction of the light scattered in polar cones has 
already been scattered by the torus. To obtain the correct net po- 
larization it is in principle necessary to conduct multi-scattering 
simulations involving both scattering regions. This is especially 
important when the base of the cone reaches far inside the funnel 
of the torus and has sufficient optical depth, and it remains true, 
even if the material in the cones is not homogeneous but orga- 
nized in localized blobs as seen in NGC 1068 (Kishimoto 19991 1. 
Furthermore, the above estimate does not include any polariza- 
tion induced by other scattering regions such as the accretion 
disk or equatorial material. 

In summary, the polarization of both type-1 and type-2 
AGNs arises in a complex way that will depend on the viewing 
angle. We predict that the polarization properties will depend on 
orientation indicators such as the radio properties, the width of 
HyS, and the presence of various types of intrinsic absorption-line 
systems. 

7.4. Polarization from the accretion disk and from equatorial 
scattering regions 

For type-1 objects, scattering from dust or electrons in the polar 
regions of the object only has a small impact on the observed 
polarization because the scattering angle is relatively small and 
both the covering factor and optical depth are low. Our model- 
ing shows that polarization induced by a dusty torus is also very 
low along type-1 viewing angles. An additional problem is that 
scattering cones only produce "type-2 polarization" (i.e., perpen- 
dicular to the symmetry axis), whilst we find that tori can only 
give the correct polarization in type-1 quasars for large open- 
ing angles. To understand the polarization properties of Seyfert-1 
galaxies, it is hence necessary to introduce additional structures 
to the common AGN scheme. 

We have reproduced the earlier result that an emission and 
scattering disk can produce type-1 polarization only when it is 
relatively optically thin regardless its geometrical thickness. The 
accretion disk itself is therefore ruled out as the source of the 
type-1 polarization. However, the disk does give intrinsic per- 
pendicular polarization that will be coupled successively with 
scattering regions such as the torus or scattering cones. 

We confirm that the correct type-1 polarization can be pro- 
duced by scattering of an equatorial disk has recently been exten- 
sively analyzed by Smith et al. ( 2002ll2004 ). The rotating BLR- 
disk they introduce together with the surrounding co-planar scat- 
tering region reproduces very well the turnover of the polariza- 
tion vector observed in the Ha and H/3 lines of many objects 
(Smith et al. l2005l l. This gives strong support for the existence of 
a flattened equatorial scattering region. Recently, additional sup- 
port for the idea comes from the discovery of a short timescale 
lag of the variations of the polarization relative to the unpolar- 
ized continuum (Shoji, Gaskell, & Goosmann 120051 Gaskell, 
Shoji, & Goosmann, in preparation). We find that the optical 
depth of the equatorial scattering material that is required to 
reproduce the observed polarization degree in type-1 objects is 
around unity, which justifies the approach of Smith et al. 2005J 

An encouraging result we obtain is that the polariza- 
tion properties at type-1 viewing angles do not depend very 
strongly on the geometry of the equatorial scattering region. 
This increases the flexibility in applying the flared disk model. 
However, according to Fig. |20]the optical depth plays a signif- 
icant role in determining the polarization. Again, it will be im- 
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portant to consider the polarization of a flared disk in relation 
to the other scattering regions and model all constituents of the 
AGN consistently. 

8. Conclusions 

We have developed a versatile Monte Carlo code for modeling 
polarization produced by scattering in many astrophysical situ- 
ations. We have shown that it reproduces previous results well. 
The code is publicly available and can be downloaded from the 
Internet. 

ModeUng polarization with a uniform-density torus for dif- 
ferent half-opening angles Oq leads to polarization degrees be- 
tween % and 20 %. The orientation of the /?- vector is per- 
pendicular to the symmetry axis for 6*0 < 53° and aligned for 
^0 > 60° at all viewing angles with i > Oq. For opening angles 
in the interval 53° < < 60°, both orientations of the polar- 
ization vector appear at type-2 viewing angles. The polarization 
behavior is not strongly sensitive to details of the torus geometry 
(e.g. cylinders with sharp edges, versus tori with smooth ellip- 
tical cross-sections). The obscuration efficiency is much higher 
for a large than for a compact torus. 

While the polarized flux of a torus scattered into edge-on 
viewing angles is rather low, we expect its polarization prop- 
erties to be relevant for the collimated light entering scattering 
cones. This is especiaUy important, if the scattering cones have 
moderate optical depths at the base. Since the opening angle of 
the torus and the cones increase with luminosity, we predict that 
the net polarization of type-2 objects depends on luminosity. 

The spectral shape of the polarization produced by scattering 
off" of a torus is nearly wavelength-independent. If a significant 
part of the the nuclear scattering in NGC 1068 is produced by 
dust in the torus then this removes the difficulty of having to 
maintain the degree of ionization of the putative electron cones. 

Polar scattering cones produce perpendicular polarization 
at all viewing angles. If they contain dust the wavelength- 
dependence of the detected radiation is different for face-on and 
edge-on viewing angles. If the dusty cones are seen face-on in 
transmission, the spectrum is reddened and the polarization is 
rather low, while the scattered spectrum at edge-on directions is 
blue and more highly polarized. 

The parallel polarization seen in type-1 AGNs arises natu- 
rally from electron scattering by optically-thin, equatorial distri- 
butions around the central source. At a given viewing angle, the 
polarization of different types of electron scattering disks does 
not depend strongly on the geometry, but it is a strong function 
of optical depth. A relatively low optical depth of only a few 
tenths can produce the observed polarization. On the other hand, 
optically thick equatorial material such an accretion disk only 
produces type-2 polarization at all viewing angles. 

All scattering regions interact, and the overall polarization 
properties depend on consistently modehng all components at 
the same time. Although the polarized flux by scattering off an 
optically-thick torus is weak at type-2 viewing angles, it has 
an important impact on the observed polarization as it colli- 
mates the radiation entering other scattering regions. The same 
is true for the accretion disk. The polarization properties of the 
individual scattering regions all strongly depend on the incli- 
nation of the system. We therefore expect, that the polariza- 
tion of AGN changes systematically with orientation parameters. 
To quantify this prediction one needs polarization modeling for 
multi-component systems. Such modeling will be discussed in a 
follow-up article (Goosmann & Gaskell, in preparation). 
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